Topological phase transitions in the non-Abelian honeycomb lattice 
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Ultracold Fermi gases trapped in honeycomb optical lattices provide an intriguing scenario where relativistic 
quantum electrodynamics can be tested. Here, we generalize this system to non-Abelian quantum electrodynam- 
ics, where massless Dirac fermions interact with effective non-Abelian gauge fields. We show how in this setup 
a variety of topological phase transitions occur, which arise due to massless fermion pair production events, as 
well as pair annihilation events of two kinds: spontaneous and strongly-interacting induced. Moreover, such 
phase transitions can be controlled and characterized in optical lattice experiments. 



PACS numbers: 37.10.Jk, 67.85.Lm, 73.43.-f, 71.10.Fd 

I. INTRODUCTION 

The ability to control systems at the quantum level offers 
an alternative and exciting route to deepen our understanding 
of nature [1]. In particular, the experimental design of many- 
body Hamiltonians provides a controlled method to study phe- 
nomena which were only believed to occur in the realm of 
condensed-matter physics. Ultracold atoms in optical lattices 
constitute a rich playground where the behavior of condensed- 
matter systems can be mimicked [2, 3], such as the superfluid- 
Mott insulator transition for bosons [4] and fermions [5, 6]. 
Often, the experimental versatility surpasses the possibilities 
of real materials, envisaging scenarios that were previously 
thought of as mere fiction. For instance, neutral bosonic or 
fermionic matter can be subjected to effective non-Abelian 
gauge fields [7, 8]. In this work, we shall describe the proper- 
ties of the non-Abelian honeycomb lattice, namely, a 2D free 
Fermi gas in a honeycomb lattice subjected to additional non- 
Abelian gauge fields. The bare Fermi gas can be considered to 
be a quantum-optical analogue of graphene [9], an interesting 
material with relativistic massless excitations that has been re- 
cently synthesized in a laboratory [10]. As shown below, the 
addition of external non-Abelian gauge fields induces a vari- 
ety of topological phase transitions beyond the usual Landau 
symmetry -breaking paradigm [11], which can be observed in 
cold-atom experiments. 

The very unusual properties of graphene, a single layer 
of carbon atoms packed in a honeycomb lattice, rely on the 
fact that the low-energy excitations display a linear disper- 
sion relation [12], and are thus described by massless Dirac 
fermions [13]. Accordingly, it is possible to observe ex- 
otic effects in low -temperature table-top experiments, which 
usually belong to high-energy physics (see [14] and refer- 
ences therein). In this context, phenomena such as Klein tun- 
neling [15, 16], or the relativistic extension of Landau lev- 
els [17-19] have already been observed. Interestingly, such 
effects might lead to the experimental realization of fully- 
relativistic Schrodinger cat states, the so-called Dirac cat 
states [20]. Besides, the transport properties of graphene are 
determined by the underlying relativistic excitations, which 



are responsible for the anomalous half-integer quantum Hall 
effect (QHE) [21-26]. These effects were also discussed 
in the context of ultracold atoms in honeycomb [9], and T$ 
(rhombic) lattices [27]. 

In this work, we describe the novel effects that take place 
in fermionic optical lattices when external gauge fields are 
switched on. Let us note that additional Abelian gauge po- 
tentials [8, 28-38] lead to analogues of well-known effects, 
such as the Hofstadter butterfly [28], the Escher staircase [29], 
the integer QHE [39, 40], the anomalous QHE in the hon- 
eycomb [41] and Kondo lattices [42], or even the fractional 
QHE [43]. When the external field is non-Abelian [7, 44, 45], 
a whole plethora of new phenomena arises. In this regime, it 
is possible to realize a quantum-optical spin Hall effect [46], 
to observe a modified Metal -Insulator transition [47], or 
to induce quasi-relativistic effects [48], such as Zitterbewe- 
gung [49]. It is also possible to observe non-Abelian ef- 
fects in atom optics [50], to obtain an optical-lattice version 
of a spin field-effect transistor [51], to ascertain the absence 
of localization in disordered relativistic one-dimensional sys- 
tems [52], or to study the effects of non-Abelian magnetic 
monopoles [53] and Yang Mills theories [54]. Such propos- 
als [7, 44] can also be used to produce effective spin-orbit 
interactions [55], to realize topological order [56, 57], or to 
generalize the integer QHE to a non-Abelian scenario [58]. 
Let us finally remark that non-Abelian features also provide 
an interesting setup where Dirac fermions emerge in a square 
lattice [40, 59], a regime previously restricted to staggered 
fields [60, 61]. 

It is the purpose of the present article to study the ef- 
fects of non-Abelian gauge fields on the emerging relativistic 
fermions that arise in the optical-lattice analogue of graphene. 
As discussed below, these fields lead to spontaneous and 
stimulated processes, where pairs of relativistic fermion- 
antifermion are created or annihilated. Intriguingly, these in- 
teracting scattering processes occur in a non-interacting Fermi 
gas, which clearly shows that the background non-Abelian 
field has dramatic consequences on the emerging low-energy 
theory. Interestingly, these creation-annihilation processes 
can be understood in the light of topological quantum phase 
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transitions [11], since the number of conical singularities and 
thus the topology of the Fermi surface is altered. Additionally, 
the creation-annihilation events occur between Dirac fermions 
with opposite topological charges [62], which allow us to fur- 
ther characterize the different topological phases. An empha- 
sis is placed upon the detection of this effect from density 
measurements. 

This paper is organized as follows, in Sec. II we review 
the properties of a 2D Fermi gas in a honeycomb optical 
lattice. In Sec. Ill, we introduce an additional non-Abelian 
gauge field, and describe its effect on the emerging relativis- 
tic fermions. Here, new phases with different Fermi surface 
topologies arise, and the gauge-induced phase transitions be- 
tween them are presented in Sec. IV. In Sec. V, we present 
a detailed study of the anomalous quantum Hall effect that 
arises in the different topological phases. In Sec. VI, we dis- 
cuss different techniques to experimentally detect and charac- 
terize these phases. Finally, Sec. VII contains our conclusions. 



II. ABELIAN HONEYCOMB LATTICE 

In this section, we review the peculiar properties of a free 
Fermi gas loaded in a honeycomb lattice (see Fig. 1(a)). In 
condensed matter, graphene provides an ideal realization of 
this system, where conduction electrons are tightly bound to 
the carbon atoms distributed along a honeycomb crystal lat- 
tice [12, 14]. The recent isolation of graphene [10] reveals an 
experimental path towards truly 2D-materials which present 
distinctive features from their 3D-counterparts. In particu- 
lar, graphene has the band structure of a semi-metal, since 
the valence and conduction bands touch at isolated points of 
the Brillouin zone. Moreover, the dispersion relation around 
such points is linear, and thus low-energy excitations behave 
as relativistic massless Dirac fermions. In accordance, this 
noteworthy material connects two different fields: relativistic 
quantum electrodynamics and condensed-matter physics. 

A quantum-optical analogue of this 2D material can be 
achieved with fermionic atoms, such as 40 K or 6 Li [5, 6], 
loaded in a honeycomb optical lattice [9] (see Fig. 1(a)). The 
corresponding optical dipole potential can be generated by 
three coplanar lasers in standing- [9] or traveling-wave con- 
figurations [63]. In fact, the bosonic superfiuid-Mott insula- 
tor transition in both honeycomb and triangular 2D lattices 
has been recently observed [64]. In order to achieve the 
half-filling and non-interacting regime, pairwise interactions 
should be controlled by means of Feschbach resonances [65]. 
Let us remark that this setup offers an optical version of free- 
standing graphene, where dislocations, impurities, phonons, 
or curved defects are completely absent. Furthermore, the 
controllability of the hopping strengths and on-site energies, 
leads to anisotropic and massive regimes that go beyond the 
real material possibilities [9, 46]. 



A. Band structure: Tight-binding model and conical 
singularities 

Let us review the band structure of fermions hopping 
along the honeycomb lattice [12] (cf. Fig. 1(a)), which 
has a two-atom basis with unit vectors ai = j(3,-\/3) and 
32 = |(3,— V^), and lattice spacing a. Fermions in the A- 
sublattice hop to three nearest neighbors of the B-sublattice 
ui = — § ( 1 , -s/3 ) , "2 = a(l,0), and 113 = §(— 1,-\/3), accord- 
ing to the following Hamiltonian 

H= -t Y,a}bj + h.c. (1) 

Here, the fermionic operators a],b\ create (annihilate) 

a fermion on i — th site of the A- and B-sublattice respectively, 
and t stands for the nearest-neighbor hopping energy. Note 
that the hopping strength is controlled through the optical po- 
tential depth [63]. The translationally invariant Hamiltonian 
in Eq. (1) can be diagonalized in momentum space by the 
introduction of sublattice operators «k = Lr J eA elkr;a .7' an d 
bk = ^ GB e lkr ^ ; -. Accordingly, the Hamiltonian becomes a 
quadratic operator H = £ k ^//k^k, where 

//k = -r£(e lku J(7 + +e- iku <c7-) (2) 

and 'I'k = (flk,^k)' is a spinor operator. Note that the Hamil- 
tonian in Eq. (2) is expressed in terms of the sublattice flip 
operators a + = | f) (I |, C7~ = ||) (f |, with t (I) representing 
the A(B)-sublattice. This Hamiltonian is readily diagonalized, 
yielding the following band structure 

E±(k) = ±t \/3 + 4 cos ¥f cos + 2 cos v^M , (3) 

which is represented in Fig. 2(a). In this figure, we observe 
that the two energy bands touch at different points inside the 
first Brillouin zone, and therefore the Fermi surface at half 
filling is reduced to a discrete set of isolated points (see also 
Fig. 1(b)). Besides, a closer inspection shows that the energy 
dispersion around these singular points is conical E ~ |k| (see 




Figure 1: (a) The honeycomb lattice can be described as two triangu- 
lar Bravais lattices A, and B, spanned by the unit vectors ai and a2- 
The nearest neighbors of the B-sublattice are A-sublattice sites given 
by ui , U2, U3. (b) The reciprocal lattice is spanned by bi , l>2, and the 
first Brillouin zone corresponds to a rotated hexagon, whose corners 
correspond to Dirac points (Note that only K, K' are independent). 
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Fig. 3(a)), which is a clear signature of the relativistic nature 
of the low-energy excitations at half filling. 

Let us now describe the content of figs. 3(a) and 3(b) in 
more detail. In Fig. 3(a), we schematically represent the dis- 
tribution of Dirac cones inside the hexagonal first Brillouin 
zone. We find six conical singularities localized on the cor- 
ners of the Brillouin zone, where the valence and conduc- 
tion bands touch (see also Fig. 2(a)). Note that only two of 
such Dirac points are independent, the remaining being ob- 
tained from this pair by mere translations in the reciprocal 
unit basis bi = §f (1, V3), b 2 = |f (1,— a/3)- The two inde- 
pendent Dirac points K, K' are represented by opaque cones, 
whereas the remaining cones are transparent. Accordingly, we 
can state that two conical singularities arise in the band struc- 
ture of the honeycomb lattice. This schematic figure is further 
supported by Fig. 3(b), where the contour of the conduction 
band in Eq. (3) has been represented. In this contour plot, 
the color scale is the following: dark blue colors represent 



a) 




Figure 2: Band structure and conical singularities of the Fermi gas in 
a honeycomb lattice, a) Graphene regime, where the bands touch 
at the corners of the hexagonal Brillouin zone. Accordingly, the 
fermion gas can be described as a semiconductor with vanishing gap 
or a semi-metal, b) ri-flux regime a = j8 = 71 / 2, where the num- 
ber of singularities is clearly raised with respect to graphene. c) 
Non-Abelian regime a = j8 = 0.2, where the valence and conduc- 
tion bands split in two, and the structure of the conical singularities 
is modified as a consequence of the coupling to the gauge field. 



energies close to zero Ef = 0, while light red colors corre- 
spond to higher energies, and thus to empty fermionic states. 
From this figure, we observe that six zero-energy points arise 
in the corners of the Brillouin zone, an hexagon represented 
with dashed lines. It is precisely at these zero-energy points 
that the conduction and valence bands touch, giving rise to the 
conical singularities formerly described. As shown in the next 
section, massless Dirac fermions emerge as the fundamental 
low -energy excitations around these singular points. 



B. Low-energy excitations: relativistic fermions 

To study the low-energy excitations, we expand the Hamil- 
tonian in Eq. (2) around the conical singularities located at 
K = |f (0, ^), K' = g (0,-^). The long-wavelength ex- 
citations are obtained after expanding the momentum around 
each Dirac point, such as K, namely p = h(k — K), which 
leads to 

Heff = ^lM P )y p ,H D =C {a xPx + OyPy) , (4) 
P 

where a x — oy, (Xy = — G x , are the Dirac matrices, which are 
reduced to the usual Pauli matrices in a 2+1 Minkowski space- 
time. This Hamiltonian is completely equivalent to the Dirac 
Hamiltonian of massless fermions, if we understand the Fermi 
velocity c = | ¥ as the effective speed of light. The rela- 
tivistic Hamiltonian in Eq. (4) readily accounts for the con- 
ical singularities E ~ c|p| observed in Fig. 2(a). Addition- 
ally, it identifies the lattice geometry (i.e. the existence of 
two interpenetrating triangular sublattices) as the responsible 
of the spinor structure of the effective relativistic theory. Let 
us note that the low energy quantum field theory associated to 
the remaining Dirac point K' also describes emerging mass- 
less fermions, but with a different representation of the Dirac 
matrices (i.e. satisfying the Clifford algebra {aj, oc^} = 2djk), 
namely, a x = G y , tt y = <3 X . Therefore, two species of indepen- 
dent Dirac fermions Nd = 2 emerge in the low-energy theory 
of a two-dimensional Fermi gas in the half-filled lattice. 

The effective relativistic picture emerging at low energies 
in Eq. (4) should be compared with the exact energy bands 
in Eq. (3). In Fig. 4(a), we represent the exact band structure 
around the Dirac point K (transparent surfaces), and compare 
it with the conical energy dispersion predicted from the Dirac 
Hamiltonian (solid surfaces). Accordingly, we show that the 
description of low-energy excitations in terms of relativistic 
massless fermions is valid in a neighborhood of the Dirac 
points (i.e. 5k = k - K < 1 /a). 



III. NON-ABELIAN HONEYCOMB LATTICE 

In the previous section, we derived the effective low-energy 
quantum field theory for the half-filled fermionic honey- 
comb lattice, and showed that two species of massless Dirac 
fermions emerge at large wavelengths, each around a differ- 
ent Dirac point. The coupling of such massless fermions to 
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Figure 3: Distribution of Dirac points with their associated topo- 
logical charges along the first Brillouin zone for (a) Graphene Pi, 
(c) .R u , -reversed graphene Pi, (e) R Ui R Ul -reversed graphene P3, (g) 
R Ul -reversed graphene P4, (i) ri-flux P5, (k) Non-Abelian point 
a = P -C 1. Contour plot of the energy bands for (b) Graphene 
Pi, (d) R U3 -reversed graphene Pi, (f) Ru 3 Rui -reversed graphene P3, 
(h) R Ul -reversed graphene P4, (i) Il-flux P5, (1) Non-Abelian point 
a = j3<l. 



external electromagnetic fields leads to interesting effects dis- 
tinctive of relativistic field theories, such as the Klein para- 
dox [16]. In this section, we shall focus on the fermion cou- 
pling to constant magnetic fields. In the Abelian scenario, 
this leads to the noteworthy anomalous integer quantum Hall 
effect, either in graphene [21-26], or honeycomb optical lat- 
tices [41]. As discussed below, this paradigm should be re- 
considered when matter is coupled to non- Abelian gauge po- 
tentials. Let us remark that even if such fields are not avail- 
able in condensed-matter experiments, it is possible to synthe- 
size them in a controlled manner by means of laser-assisted 
tunneling processes in optical lattices [7], or dark-state meth- 
ods [44]. Therefore, surpassing graphene-based materials, op- 
tical lattices present exotic avenues at the forefront of con- 
densed matter, high-energy physics, and atomic physics. 

In non-Abelian lattice gauge theory, a unitary matrix Un 
is associated to the link connecting the lattice points r, — > ry, 
which must be generated by the elements of the underlying 
Lie group. Consequently, the Abelian hopping Hamiltonian 
in Eq. (1) must be generalized to 



H 



(hJ) 



Ultimo 



h.c, 



(5) 



where a' iT , bj T (a iz , bi x ) create (annihilate) a fermion at lat- 
tice site r, of the A- and B-sublattice, with a color-like com- 
ponent T = 1,2, ...,N corresponding to the fundamental rep- 
resentation of the Lie group. Such a situation has been de- 
picted in Fig. 5 for the particular case of SU(2) fields (i.e. 
fermions have two colors T = 1,2 ). Here, we observe that 
the coupling of matter to non-Abelian fields induces unitary 
transformations along each hopping path: U\ = e lttT> , U2 = I, 
and C/3 = e^ Tv , where T x ,T y are Pauli matrices in the color- 
components, and a,j3 represent the non-Abelian fluxes (see 
Fig. 5). In the momentum representation, the Hamiltonian in 
Eq. (5) becomes H = £ k M^Z/k^k, where 



e iku ^ 



)^+e- iku ^- ( 



(6) 



and the spinor ^Pk = («ki , flk2, £>ki , bk2)' includes pseudo-spin 
and color degrees of freedom, which correspond to the under- 
lying sublattices and to the internal fermionic components. 
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Figure 4: Comparison between exact energy bands and the low- 
energy approximation, a) In graphene, we represent the exact en- 
ergies in Eq. (3) (transparent colors) and the energy spectrum from 
the low-energy Hamiltonian in Eq. (4) (solid colors) around K. b) 
Idem for a = /3 = 0.2, comparing Eq. (25) to Eq. (4) around K; . 
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In the following sections, we study the pattern of emerg- 
ing relativistic fermions associated to different gauge fluxes 
(a, j3). As shown below, the total number of conical singular- 
ities and their distribution within the Brillouin zone depends 
on the particular values of the gauge fluxes. We shall char- 
acterize different Abelian and non-Abelian phases in terms of 
these low-energy massless fermions. 

A. Wilson loop and non-Abelian regimes 

The non-Abelian gauge fields that induce these conditional 
hoppings are defined through a generalized Peierls substitu- 
tion Uij = V^e'^^^Vj, where the independent matrices Vj 
reflect the local gauge invariance, and the gauge field A^(ri) 
is located at the link ri = j(r,- +rj) and directed towards 
fl = Yj — r,. The particular choice U\ = e iaT \ C/ 2 = I, and 
E/3 = e'^ Tv in Fig. 5, is locally equivalent to the following 
SU(2)-generated gauge field A = La^^^( r i) T « e /j 

A = a(T x ~ Xy)e x + ~ (t x + x y )e y . (7) 

In the continuum limit a — > 0, and according to non-Abelian 
gauge theory, the associated strength tensor F^ v = d^Ay — 
<5vAfi — ifA^Av], describes a constant magnetic field along 
the z— axis B = ^=-T z e z . On the lattice, these external fields 
are further characterized by the Wilson loop W = trf/o, where 
Uq = U1U2U3U1U2U* is the unitary transformation around 
an elementary plaquette. The Wilson loop represents a gauge- 
invariant non-local observable related to the non-Abelian flux 
through an hexagonal plaquette, and is expressed as follows 

W(a,P) = 2cos 2 a + 2cos2j3sin 2 a. (8) 

In Fig. 6, we depict the Wilson loop as a function of the exter- 
nal fluxes, and identify the Abelian (|W| = 2) and non-Abelian 
(\W\ 7^ 2) regimes. 




Figure 5: Scheme for the fermionic honeycomb lattice subjected to 
SU(2) gauge fields, where each hopping is dressed by Ui = e iaTj: , 
U2 = I, and I/3 = e'P^'. We have used double-line links to indicate 
the underlying color degrees of freedom % = 1 , 2 associated to each 
fermion. The gauge fields on the links induce a color-dependent hop- 
ping between adjacent sites, which has been represented by squares. 



Abelian regimes: The following values of the gauge fluxes 
correspond to Abelian phases in the Hamiltonian of Eq. (5) 

• Graphene point Pi for (a,j3) = (0,0), 

• R U3 -reversed graphene P2 for (a,j8) = (0,n), 

• R U3 R Ul -reversed graphene P3 for (o£,/J) = {it,n), 

• R Ul -reversed graphene P4 for (a,j3) = (ti,0). 

• n-flux point P 5 for (a,j3) = {\n, \n). 

Note also that the lines connecting Pi,Pi,P$,P4 correspond 
to Abelian phases (see Fig. 6). In the next section, we de- 
rive the band structure and low-energy massless fermions for 
these different regimes. We show that the four graphene-like 
points P\ , P2, P3 , Pa have Nj = 4 species of emerging Dirac 
fermions (notice the doubling due to color-degeneracy), and 
are thus topologically equivalent. Conversely, the Il-flux 
phase presents N e i = 8 different gapless fermions, and thus a 
distinct Fermi surface topology. As shown in Sec. V, transport 
properties in this Il-flux phase, such as Hall conductances, 
clearly deviate from the graphene-like phases. 

Non-Abelian regimes: The remaining values of the gauge 
fluxes define non-Abelian phases with \W\ < 2, whose low- 
energy properties shall dramatically depend on the particular 
values of the fluxes. Below, we study the non-Abelian regime 
a = j3 <C 1 close to the Graphene point P\, where the low 
energy dispersion is linear and saddle-point in orthogonal di- 
rections. In this respect, the effective mass of quasiparticles 
at large wavelengths is highly anisotropic, and leads to novel 
features in the transport properties. 

B. Effective low energy theories in Abelian regimes 

In this section, we study the emerging quantum field theo- 
ries of the non-Abelian honeycomb lattice at half-filling, cen- 
tering our attention to the Abelian regimes introduced above. 

Graphene point: In the limit of vanishing gauge fluxes 
a = = (see point Pi in Fig. 6), we recover the standard 



i'3 




Figure 6: Wilson loop associated to the hexagonal plaquette under 
the gauge fields in Eq. (7). The Abelian regimes (\W\ = 2), corre- 
spond to the flux configurations P\ = (0, 0), Pi = (0, K), F3 = (it, it), 
P\ = (?r,0), the straight lines joining them, and the marginally 
Abelian point P5 = (5, j). Different paths, such as Yi and Y2> j°i n 
the above regimes by means of a quench in the external fluxes. 
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graphene Hamiltonian of Eq. (2), with an additional 2-fold 
color degeneracy 

H k = -t (e iku ' + e lkU2 +e ikU3 ) ff+ ®I 2 +h.c. (9) 

Therefore, energy bands correspond to those expressed in 
Eq. (3) and shown in Fig. 2, with the peculiarity that each 
band is doubly-degenerate. As a consequence, the number of 
independent Dirac points is raised to Nd = 4, fulfilling 

Hg = c(Oy®hp x -a x ®hpy), (10) 
around the singularity K, and 

H§' = c(Oy® hp x + Ox®*2Py) , (11) 

around K'. As expected, the species of emerging massless 
fermions is doubled due to the color degrees of freedom. 

R„ 3 -reversed graphene: In the Abelian phase correspond- 
ing to a = 0, [5 = n (point P2 in Fig. 6), the Hamiltonian 

H k = -t (e lkui + e ikU2 -e ikU3 ) a+ ®I 2 +h.c, (12) 

can be obtained from graphene's analogue in Eq. (9), after 
reversing the hopping amplitude along 113 (i.e. t — > —t). The 
energy bands fulfill 

E± (k) = ±t \fi - 4 sin sin - 2 cos V3k y a , (13) 

whose contour clearly shows that the pair of conical singu- 
larities are no longer at the corners of the Brillouin zone, but 
well inside it (see the contour energy in Fig. 3(d)). The pair 
of conical singularities, now located at K = 2^75)' an( l 

= If ("ST' 273^' correspond t0 tne corner Dirac points of 
graphene, which have been translated in momentum space 
along 120° direction. Such interesting transport of Dirac 
cones around the Brillouin zone is schematically represented 
in Fig. 3(c). Here, transparent cones represent the initial Dirac 
fermions for vanishing fluxes a — j3 = 0, and opaque cones 
represent the transported cones at a — 0, j3 = %. The direc- 
tion of transport is shown by black arrows. The effective low 
energy theory at half-filling around these points is also cap- 
tured by the relativistic Dirac equation, with color-degenerate 
massless fermions (i.e. a total of Nj = 4 fermions), that be- 
have according to the Hamiltonian 

H% = c{a x ®l 2 p x + ay®l2Py), (14) 

around the singularity K, where a different representation of 
the Clifford algebra is chosen a x = —& y ,(Xy = G x , where 

/ w\ _ / -iw\ 
°* = { w* ) ' a >={i W * ) ' (15) 

and w = 2(1 — i\/3) lies in the complex unit circle. Around 
K', the low-energy Hamiltonian leads to (X x = 6 y ,(Xy = O x . 



R U3 Rui -reversed graphene: Considering the gauge fluxes 
a = [5 = K corresponding to the point P3 in Fig. 6, the Hamil- 
tonian 

H k = -t (-e lkui + e ikU2 -e ikU3 ) a+ ®I 2 +h.c, (16) 

is analogous to graphene with the reversed hoppings t — > — t 
along the directions ui , U3. The underlying energy bands are 

£±(k) = ±ry / 3-4cos^'cos^|^+2cosv / 3M, (17) 

and lead to the conical singularities depicted in figs. 3(e) and 
3(f). The pair of conical singularities lie at K = |f (0, A^), 

and K' = §f (0, 5|), an d correspond to the corner singularities 
translated in momentum space along a 180° direction. At half- 
filling, the low energy excitations around such points are also 
relativistic massless fermions 

Hg = c(a x ® I 2 p x + Oy ® hp y ) , (18) 

with a x = G y , (Xy = a x for K, and a x = o y , a y = —a x for K'. 

R u , -reversed graphene: Considering the gauge fluxes a = 
7T,/3 = corresponding to P4 in Fig. 6, the Hamiltonian 

H k = -t (-e lkui + e ikU2 +e ikU3 ) a + ®I 2 +h.c, (19) 

is analogous to graphene with the reversed hoppings t — > — t 
along the direction Ui. The underlying energy bands become 

E±(k) = ±t \/3 +4 sin ¥f sin - 2 cos y/3k y a, (20) 

represented in figs. 3(g) and 3(h). The pair of singularities lie 
at K = f a {=±, ^), and K' = f a {\, ^|), and correspond to 
the corner singularities translated in momentum space along 
a 60° direction. At half-filling, the low energy excitations 
around such points are also relativistic massless fermions 

Hn=c(a x ®l2p x + a y ®Izp y ), (21) 

with a x = —Oy, and a y = O x for K, where &j correspond to 
Eq. (15) for w = i(l +iy3). Equivalently, the K' emerging 
fermion leads to the Clifford algebra a x — & y , and a y = a x . 

Let us finally remark the fact that all the Abelian phases 
described so far simply correspond to finite translations of 
the graphene points, and thus present the same number of 
emerging fermions N c / = 4. Accordingly, the physical prop- 
erties of the two-dimensional Fermi gas at half filling are ex- 
pected to be analogous to graphene. The unique effect of a 
non-vanishing Abelian flux a,j5 is to transport the relativis- 
tic fermions around the Brillouin zone. Below, we study 
richer situations where the topology of the Fermi surface is 
modified by a variation of the number of gapless fermionic 
species. It is also important to note that these Abelian phases 
can also be achieved by simply tuning the anisotropy of the 
hopping parameters [9], without the need to insert external 
fields. Nonetheless, the phases to be discussed below cannot 
be achieved in this fashion, and constitute a non-trivial conse- 
quence of the applied gauge fields. 



7 



Il-flux point: Let us now focus on the Abelian regime a = 
j3 = | , which corresponds to the II-flux phase (i.e. lattice 
fermions pick up an overall minus sign by hopping around the 
elementary plaquette). This regime corresponds to the point 
Ps in Fig. 6, where the fermionic Hamiltonian 



H k = -t (^ie ikui cr + ® x x + e 1KU2 a + ® I 2 + ie 1KU3 a + ® tJ + h.c. 

(22) 

induces a color flip transformation when the fermion hops 
along directions Ui and 113. The energy bands 



Jku 2 



iku 3 rs+ , 



£(k) = ±fy3±y6 — 4cos3£ T acos \f%k y — 2cos2y/3k y a, 

(23) 

have been represented in Fig. 2(b), where one observes how 
the number of conical singularities is increased with respect 
to the previous Abelian phases. Let us note that the color- 
degeneracy is lifted due to the color-flip induced hopping. The 
distribution of conical singularities in the Brillouin zone can 
be better appreciated in figs. 3(i) and 3(j), where N c [ = 8 inde- 
pendent Dirac points appear in positions 
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(24) 



We conclude that the topology of the Fermi surface in the 
Il-fiux phase is not equivalent to the remaining four Abelian 
phases described so far. Indeed, the number of Dirac fermions 
is distinct, which reveals the existence of processes, deeply 
rooted in the non-Abelian phase, where fermions scatter and 
pair-production (annihilation) is allowed. Such processes 
shall be discussed in detail in the forthcoming section, where a 
topological argument will clarify why such interacting events 
occur even when the Fermi gas is non-interacting. 



C. Effective low energy theories in a non-Abelian regime 

In order to identify the effects of the non-Abelian fields, we 
focus in the regime a = j5 -C 1, where the topology of the 
Fermi surface can be worked out analytically. We note that 
a variety of different and interesting non-Abelian phases also 
occur for different fluxes ( a , j5 ) , leaving their numerical treat- 
ment for the following sections. The energy band structure 



E (k) = ±t y f a ,k ± s/gafr, 
can be expressed in terms of the following functions 



(25) 



fak =3+2or +4cos 
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-2cosv3&y) 



ga.k = - 2a 2 (2 cos 3^, cos \/3k y + (2 + a 2 )cos2\/3k y (26) 
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These energy bands have been represented in Fig. 2(c) for the 
particular case of a = 0.2, where one readily observes that 
the non-Abelian features of the external fields dramatically 
modify the corresponding conical singularities. The valence 
and conduction bands each split into two components, and the 
conical singularities are displaced from the corners of the Bril- 
louin zone (see also Fig. 3(k) and 3(1)). However, we note that 
the two initial (two-fold degenerate) Dirac points give birth to 
four (non-degenerate) independent Dirac cones, so that Nd = 4 
remains constant in this non-Abelian regime. The exact posi- 
tion of the displaced Dirac points depends on the non-Abelian 
flux a as dictated by the following expressions 
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(27) 



which reduce to the Abelian points for a = as expected. 
Once the exact location of the non-Abelian Dirac points has 
been identified, it is possible to linearize the Hamiltonian 
in Eq. (6) and obtain the continuum theory describing low- 
energy excitations. In the particular case of p = h8k = 
h(k. — Ki), we get the free massless Dirac Hamiltonian 



H£fac(a x ®hPx 



l 2Py )+AH, 



(28) 



with a perturbation AH describing spin- and color-flip pro- 
cesses that arise due to the non-Abelian fluxes 

AH = Ai a + ®I + A 2 CJ + (g)T + +A 3 cj + (g)T _ +h.c. (29) 

Here, we have also defined 

^ = (1+2$) - (V3(l-§ 2 )- 4 )5*y + i(l -t;)$k x , 
A 2 = aV2ri + - 

A3 = —aV2rj- 



^^Sky+^Sk*, 



^Brj+Sky + ^rj-Skx, 



(30) 



with the following parameters rj± = ± y/ 1 — $ 2 )e ljr//4 , and 
§ = (-1 + Va 2 + 2a 4 )/2(l + a 2 ). From the expressions 
above, one can check that in the vanishing flux limit a — > 0, 
the perturbation also vanishes AH —> 0, and we recover the 
Abelian graphene theory. Let us remark here that the non- 
Abelian perturbation is essential to account for the displace- 
ment of the Dirac points. As follows from Fig. 4(b), the dis- 
persion relation around the singular points is not simply con- 
ical. Indeed, at low energies, the spectrum becomes linear in 
one direction, and displays a saddle point in the transverse di- 
rection. In this respect, we can expect a dramatic modification 
of the massless fermions properties due to their coupling to 
these exotic non-Abelian fields. 
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We conclude that the effect of the synthetic gauge fluxes is 
two-fold: they modify the Fermi surface by transporting the 
Dirac points according to Eq. (27), and they are also respon- 
sible of a different low-energy theory in Eq. (28). Indeed, 
the specific form of the perturbation in Eq. (29) suggests the 
existence of certain emerging gauge fields that couple to the 
bare massless fermions and give rise to this two-fold effect. A 
similar effect was explicitly shown in graphene-like structures 
with a Kekule distortion [66, 67]. 

IV. TOPOLOGICAL QUANTUM PHASE TRANSITIONS 

In the previous sections, we have thoroughly described the 
pattern of Dirac points and emerging relativistic fermions for 
different Abelian and non-Abelian phases. These phases dif- 
fer by the total number of Dirac fermions and their distribu- 
tion within the Brillouin zone. Below, we characterize the 
topological properties associated to the pattern of massless 
fermions by means of a topological charge (i.e. topological 
invariant) ji = ±1 defined for each conical singularity [62]. 
Moreover, the topological properties of the Fermi surface for 
non-interacting Fermi systems fully determine the underly- 
ing quantum order, the paradigm of order in quantum states 
beyond the Landau symmetry -breaking description [11, 68]. 
In this respect, the pattern of topological charges classifies 
the phases introduced above within different quantum-ordered 
universality classes having the same symmetry. 

We also study how a modification the gauge fluxes (cc,j3) 
connects different quantum orders (i.e. different patterns of 
topological charges). In this respect, non-Abelian gauge fields 
induce a quantum phase transition (QPT) (i.e. a phase transi- 
tion at T = driven by quantum fluctuations) between differ- 
ent vacua with the same symmetry, and cannot be classified by 
the Landau symmetry-breaking paradigm. This class of QPTs 
have been previously studied in gapped systems such as the 
fractional QHE [69], or gapless strongly correlated systems, 
such as certain types of spin liquids [70, 71]. In the sections 
below, we show how such QPT occurs in the realm of Fermi 
gases subjected to non-Abelian fields, which can be probed in 
a cold atom experiment. 

A. Topological characterization: winding numbers 

The topology of the Fermi surface determines the type of 
effective theory and emergent symmetries at low energies. In 
order to classify the different universality classes in Fermi sys- 
tems with gapless excitations, a topological winding number 
can be defined for the general Hamiltonian H = J^. ^//k^Pk, 
where was defined in Eq. (6). This Hamiltonian fulfills a 
particle-hole symmetry {F,^} = 0, where F — diagQk, — 12), 
and allows us to define the following winding number 

V = 4^£ tr ( ra k 1<ZH J (3D 

where C is a curve defined in momentum space where the 
Hamiltonian presents no zeros [62]. Such a winding num- 



ber, or the equivalent version 

ii = — I <il(det/z k ) _1 <9idet/ik, (32) 
2ni Jc 

where = — t L ; -e~ lku Jf/J, has three important properties: 

Dirac fermions: The zeros of H^, and thus conical singu- 
larities, are detected inside the region delimited by the curve 
C, by means of a non-vanishing winding number jj, 7^ 0. 

Homotopic invariant: The winding number is a homo- 
topic invariant, and is thus independent with respect to dis- 
tortions of the loop C as far as H^ 1 exists. 

Topological invariant The winding number is a topolog- 
ical invariant, and is thus independent with respect to local 
perturbations of the Hamiltonian Hk — > + SH^. 

At half-filling, the Fermi surface is reduced to a set of 
points and we can characterize its topology by the pattern 
of conical singularities and their corresponding topological 
charges (K;, /J/). As discussed in the next section, varying the 
gauge fluxes (ct,j3) along a non-Abelian path 7: (ao,j3o) — > 
(ttf,[lf), such that \W\y ^ 2, we observe that different quan- 
tum orders are connected by quantum phase transitions. We 
have first determined the pattern of topological charges in each 
phase, and the corresponding universality classes: 

Graphene universality class: The quantum order in this 
equivalence class is determined by two positive topological 
charges ji = +l, and two negative topological charges ji = — 1 
within the first Brillouin zone (see Fig. 3(a)). There are thus 
N c i = 4 Dirac points, which can be distributed in different po- 
sitions in momentum space (see figs. 3(a)-3(h)), maintaining 
the same quantum numbers and thus belonging to the same 
universality class. We can thus conclude that all the graphene 
points Pi,P2,Pi,P4 describe the same quantum order. 

Li-flux universality class: The quantum order in this 
equivalence class is determined by four positive topologi- 
cal charges ji = +1, and four negative topological charges 
/i = — 1 within the first Brillouin zone (see Fig. 3(i)). There- 
fore, a total of N c [ = 8 species of gapless fermions arise. 

As discussed so far, the different quantum orders are all de- 
termined by a distinct even number of gapless fermions. The 
adiabatic variation of the external gauge fluxes will give rise 
to processes, where the number of massless fermions is mod- 
ified. Note however, that the total number of massless exci- 
tations shall always remain even, a distinctive feature of all 
the phases discussed above. A different prescription occurs 
for topological insulators with time-reversal invariance, where 
the number of massless excitations can be generally odd [72- 
81]. We believe that it would be interesting to study this type 
of topological phase transitions in such systems. 

B. Abelian path: absence of a quantum phase transition 

In this section, we study the Fermi surface topology as the 
gauge fluxes are modified along an Abelian path 7: (ao, j3o) — > 
(af,Pf) such that \W\y = 2. According to the previous sec- 
tion, this path only joins regimes belonging to the same uni- 
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versality class, that of graphene. Therefore, no phase transi- 
tion occurs, and modifying the gauge flux only transports the 
Nd = 4 topological charges in momentum space. We focus on 
the following Abelian path 71 (f ) = (tn,0), with t G [0, 1], con- 
necting the 1 st graphene point Pi to the 4 th graphene point P4, 
as represented in Fig. 6. In figs. 7(a)-7(h), the distribution of 
conical singularities along the Abelian path 71 is represented. 



a) b) 




Figure 7: Pattern of Dirac points along the Abelian path connecting 
graphene P\ to 4 th graphene point P4. The gauge fluxes are the fol- 
lowing (a) a = 0, j8 = 0, (b) a = 0. 15s, /3 = 0, (c)a = 0.30s, j3 = 0, 
(d) a = 0.45s, /3 = 0, (e) a = 0.60s, /3 = 0, (f) a = 0.75s, j8 = 0, 
(g) a = 0.90s,/3 = 0, (h) a = s,j3 = 0. 

In these figs. 7(a)-7(h), the Brillouin zone has been de- 
picted by a dashed hexagon, and the conical singularities cor- 
respond to the dark-blue points where the valence and conduc- 
tion bands touch. This set of contour figures represent the con- 
duction band energy obtained by diagonalizing the Hamilto- 
nian in Eq. (5) for different fluxes along the path 71 . It clearly 
follows that, starting from graphene Pi in Fig. 7(a), the initial 
doubly-degenerate topological charges are split into a pair of 



Dirac points that travel in opposite directions (see Fig. 7(b)). 
Note however that the total number of Dirac points Nj=A and 
topological charges ji = {— 1,— 1,+1,+1} are always pre- 
served from Fig. 7(a) to Fig. 7(h). Having the same pattern of 
topological charges, the different regimes along the Abelian 
path all belong to the same universality class, and have the 
same quantum order. The effect of a varying gauge flux along 
this path is simply the transport of topological charges in mo- 
mentum space. Therefore, no QPT connecting two different 
quantum orders exists, a fact that also holds for any other 
Abelian path (e.g. that joining P2 to P3 in Fig. 6). 



C. Non-Abelian path: topological quantum phase transition 

From the previous section we conclude that the Abelian 
regime must be abandoned if different quantum orders are 
to be connected. Below, we show that a non-Abelian flux 
quench provides a rich scenario where scattering processes 
occur, involving the creation or the annihilation of massless 
Dirac fermions. Moreover, the pattern of topological charges 
is modified accordingly. Let us emphasize the exotic nature of 
these interaction processes, which occur in a non-interacting 
gas. As shown below, only scattering of fermions with oppo- 
site topological charges is allowed. Moreover, we can draw 
an interesting analogy for two types of topological QPTs, 
which can be driven by a stimulated or spontaneous creation- 
annihilation of opposite topological charges. 

Stimulated creation-annihilation of Dirac fermions: 

Here, we study a non-Abelian path 72 that joins the 1 st 
graphene point Pi to the II-flux regime P5, and then to the 
3 rd graphene point P3 (i.e. 72(f) = (t7C,t7C), with t e [0, 1], as 
represented in Fig. 6). This path is clearly non-Abelian, since 
\W\yj =/=2 except at t = \ where the II-flux regime is reached, 
and, as shown below, the pattern of topological charges deter- 
mining the underlying universality class is abruptly modified 
for certain critical fluxes (a c ,/3 c ). 

In figs. 8(a)-8(h), we observe how the initial doubly- 
degenerate Dirac points (Fig. 8(a)) are split into a pair of 
fermions which move along the vertical axis in opposite direc- 
tions (Fig. 8(b)). Remarkably, in Fig. 8(c) we can observe how 
additional Dirac fermions are being created along the borders 
of the Brillouin zone. A closer inspection of the energy band 
dispersion around k = §f (— 1 , 77=) (i.e. the region where these 
gapless fermions are created), reveals that there is a process 
of induced or stimulated pair production (see figs. 9(a)-9(d)), 
where a couple of massless fermions is created in the pres- 
ence of a third gapless mode. Interestingly enough, the pair of 
created fermions carry opposite topological charges \i = ±1, 
so that the total winding number is conserved during the scat- 
tering process (see also Fig. 9(e)). Note that the production 
of fermionic gapless modes occurs at T = 0, and thus a QPT 
between different quantum orders takes place. In the same 
manner symmetry-breaking phase transitions lead to gapless 
Goldstone bosons, QPTs driven by fluctuations of quantum 
order can give rise to gapless fermionic excitations. Also, 
analogous to the symmetry protection of Goldstone modes, 
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these massless fermions are protected by quantum order (i.e. 
topological invariant charges jj. = ±1). 

The quantum phase transition taking place around a c = 
[5 C = |, connects two different phases with Nj = 4 and Nd = 8 
low energy Dirac points, so that the total topological charge 
is conserved. Let us also note that there is a critical point 
at a L = f$ c = x wnere tne inverse process takes place. In 
this case, two massless fermions with opposite charges are 
annihilated and connect the phases with Nd = 8 and Nd = 4 
fermions. It is also interesting to remark that at criticality 
<X C = P c — f > one can show that the massless excitations are 
highly anisotropic, being relativistic in one direction, and non- 
relativistic along the transverse axis. 

Spontaneous creation-annihilation of Dirac fermions: 

Below, we prove that spontaneous processes of pair creation- 
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Figure 8: Pattern of Dirac points along the Abelian path connecting 
graphene P\ to the 3 ld graphene point P$, passing through the IT-flux 
regime P5. The gauge fluxes are the following (a) a = = j8, (b) 
a = 0.15s = j3, (c)ce = 0.30s = j8, (d) a = 0.50s = j8, (e) a = 
0.60s = p, (f) a = 0.75s = j3, (g) a = 0.90s = j3, (h) a = s = J3. 



annihilation are also possible, where the pair of oppositely 
charged fermions are created from a local vacuum (i.e. locally 
gapped spectrum). Here, we study a non- Abelian path 73 that 
joins the 1 st graphene point P\ to the (a, j3) = (s/2,s/4) (i.e. 
Ji{t) = (t7C,tn/2), with t € [0,1]). The different Dirac points 
along such a path have been represented in figs. 10(a)-10(h). 
In Fig. 10(d), we observe how a pair of fermion-antifermion is 
spontaneously created from a locally gapped vacuum around 
two corners of the Brillouin zone k = gf f° r a — 

2/3 = 0.6s. Then in Fig. 10(e), we observe how one of the re- 
cently created fermions scatters with another Dirac point and 
how they are spontaneously annihilated for a — 2/3 = 0.7s. 

In figs. 1 l(a)-l 1(d) we study the energy band landscape as 
the gauge fluxes are swept across the critical point a c = 2/3 6 = 
0.63s. Initially, the low energy excitations are clearly gapped 
and for low enough energies we have a local vacuum. Then, 
as the gauge fluxes increase, we observe how the local gap 
vanishes and two massless fermions appear. The topologi- 
cal charge distribution is schematically shown in Fig. 11(d), 
where the total topological charge is once more conserved. 
Let us note that for a c = 2/3 c - = 0.7s (see Fig. 1 1(f)), the time- 
reversed scattering process takes place. Here, we observe a 
pair fermion-antifermion being annihilated due to their oppo- 
site charges. These events of spontaneous pair creation or an- 
nihilation connect two distinct quantum orders Nd = 4 — > 8 
by means of a non-symmetry breaking quantum phase tran- 
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Figure 9: Energy bands dispersion around the k-space region where 
an induced pair production takes place. The gauge fluxes are the fol- 
lowing (a) a = J3 = 0.2s, (b) a = J3 = 0.25s, (c)ce = /3 = 0.3s, (d) 
a = j3 = 0.35s. ( e) Scheme of the process of induced pair produc- 
tion representing the phenomena in Fig. 8(c), and the reverse annihi- 
lation process (time-reversed event) corresponding to Fig. 8(f). 
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sition. A similar spontaneous creation-annihilation process is 
also illustrated in fig. 16 (see Sec. V). In the following sec- 
tions, we address the issue of characterizing the above topo- 
logical phases through accessible measurement techniques in 
the optical-lattice setup. 

Let us close this section by comparing the topological QPT 
presented in this work, to the semi-metal to insulator transi- 
tion in bare graphene-like systems [9, 82-87]. In graphene, a 
modulation of the hopping anisotropy induced by strain can 
lead to the merging of Dirac points, which turns graphene 
into an insulator. An insulating gap can also be dynamically 
induced by quenching the sublattice energy imbalance [86]. 
Conversely, in the non-Abelian honeycomb lattice, the merg- 
ing of Dirac points does not lead to a global gap, and the semi- 
metallic phase is always preserved. Accordingly, this optical- 
lattice setup offers a broader spectrum of topological QPTs, 
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where both spontaneous and stimulated creation-annihilation 
of Dirac fermions occur. 



V. THE QUANTUM HALL MEASUREMENTS 

The presence of several distinct Dirac points within the first 
Brillouin zone leads to dramatic consequences in the physi- 
cal properties of the system. In the low-energy regime, these 
relativistic effects are adequately described by the underly- 
ing effective Dirac theories. An interesting effect to be con- 
sidered in this context is the so-called anomalous quantum 
Hall effect (AQHE) [21, 22], which has been discovered in 
graphene sheets subjected to external magnetic fields [23- 
26]. This effect is characterized by the peculiar sequence 
of plateaus in the transverse conductivity of the system as a 
function of the chemical potential: in the low-energy regime, 
where the spectrum exhibits the so-called relativistic Landau 
levels [17-19], the Hall conductivity is quantized according to 
Oxy = ±g (V + I), where v is an integer, h is Planck's constant 
and we set e = 1 . The Hall measurements show steps of four 
that can be traced back to the presence of two independent 
two-fold degenerate Dirac points N,/ = 4, and a half-integer 
anomaly due to the zero-energy mode associated to the self- 
conjugate lowest Landau level. 

In previous sections, it has been shown how several Dirac 
points are induced in our cold-atom system. In particular, we 
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Figure 10: Pattern of Dirac points along the non-Abelian path con- 
necting graphene P\ to a = 7T,j8 = j7t. The gauge fluxes are the fol- 
lowing p = \a (a) a = 0, (b) a = 0A5n, (c)a = 0.3rc, (d)a = 0.57t, 
(e) a = 0.67T, (f) a = 0.7tt, (g) a = O.Sit, (h) a = n. 



Figure 1 1 : Energy bands dispersion around the k-space region where 
a spontaneous pair production takes place. The gauge fluxes are the 
following (a) a = 2j3 = 0.67T, (b) a = 2j3 = 0.63?r, (c)ce = J3 = 
0.647T, (d) a = 2/3 = 0.667T. (e) Scheme of the process of sponta- 
neous pair production, and annihilation (time-reversed event). 
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have shown that the number of independent Dirac cones N c / 
is not necessarily conserved as the non-Abelian fluxes ((X,j5) 
are varied. The aim of this section is to show that the number 
of emerging massless fermions N c / can be inferred from Hall 
measurements. In particular, we show that the Hall plateaus 
at zero Fermi energy Ep = fulfill 

2h 

and thus provide a Dirac fermion witness. Moreover, the 
robustness of this witness is associated to the system zero- 
modes, and is guaranteed by the Hamiltonian particle-hole 
symmetry and topological index theorem [21]. The relation 
of the Hall conductivity to the number of Dirac fermions in 
Eq. (33) can be conjectured from the theoretical results in 
[21-24]. Below, we confirm this conjecture by numerically 
computing the Chern numbers associated to the energy band 
structure, and thus estimating the corresponding Hall conduc- 
tivities. Furthermore, this conductivity can be experimentally 
accessed in cold atoms by the so called Streda formula [88]. 

In this section, we shall first review how such measure- 
ments can be performed in cold-atom experiments by analyz- 
ing the density profiles. Then, we show how the AQHE taking 
place in this system gives us crucial information on the topo- 
logical QPTs that occur as the non-Abelian fluxes are varied. 
In particular, we discuss under what conditions the Hall mea- 
surements are indeed witnesses of such transitions. 



A. The anomalous quantum Hall effect and the density profiles 

In order to provoke the quantum Hall effect in our system, 
one has to open the central gap. This can be achieved by sub- 
jecting the system to an additional Abelian gauge field (i.e. an 
artificial magnetic field). Such synthetic magnetic fields can 
be produced for neutral atoms by rotating the lattice [33-37], 
by considering laser-induced methods [8, 28-32], or even by 
immersing the system in a rotating BEC [38]. The presence 
of an external Abelian gauge field modifies the hopping oper- 
ators introduced in Sect. Ill, which now take the form 

Vi =e a*4* e ta<% U 2 = I, U 3 =e i ^ (34) 

where <f> is the effective magnetic flux quanta per unit cell. A 
clear manifestation of the external magnetic field can already 
be found by computing the energy spectrum: when the flux is 
given by the ratio of two integers, namely <t> = p/q, the initial 
band structure is split into several subbands. For such rational 
Abelian fluxes, we note that the Hamiltonian takes the form 
of a Aq x 4q matrix. Therefore the number of gaps highly 
depends on the integer q, but also on the specific values of 
the non-Abelian fluxes a and j3 . When plotted as a function 
of the flux 4>, these numerous subbands form a fractal 
structure which reminds the well-known Hofstadter butter- 
fly [89], as depicted in Fig. 12 for the Abelian case a = j3 = 0. 
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Figure 12: (top) The Hofstadter butterfly for an hexagonal lattice 
subjected to a magnetic flux <5. Here one sets a = jS = 0. (bottom) 
The same butterfly spectrum in the vicinity of the Dirac regime E £ 
[— 1, 1]. The main energy gaps at <J>i = 1/157 and $2 = 2/157 are 
labelled with different symbols. 

In the cold-atom framework, the conductivity tensor mea- 
sures the response of the system to an external perturbation, 
such as a lattice acceleration, which generates a current. In 
the presence of an additional magnetic field, a Hall current is 
generated and this transverse transport is dictated by the Hall 
conductivity o#. When the Fermi energy Ep lies in a spec- 
tral gap, the Hall conductivity can be expressed as a sum of 
topological invariants 

°&=fc L N ch(En), (35) 
n E„<E F 

where N c h(E n ) is the Chern number associated to the filled 
energy band E n , which is necessarily an integer [90]. Con- 
sequently, the sum in Eq. (35) characterizes the different 
topological phases associated to the spectral gaps. These 
quantum Hall phases can be reached by varying the Fermi 
energy Ef, namely, by controlling the atomic filling factor. 
The computation of the Chern numbers can be achieved 
numerically through diverse techniques [40, 58, 91] and 
leads to the complete description of the QHE in our system. 



13 




Figure 13: (a). Hall plateaus in the low-flux regime $<C 1 and 
a = P = . (b) The density of states in the same configuration. 



In Fig. 13 (a), we show the Hall plateaus in the Abelian 
case a = j3 = 0, as recently discovered in graphene. One 
clearly distinguishes the Dirac regime, situated in the range 
Ep £ [—1,1], where the anomalous double steps and a 
half-integer plateau are observed. It should be stressed that 
the range of this relativistic regime, where the dispersion 
relation is linear, is known to be bounded by singularities in 
the density of states (see Fig. 13 (b)). 
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Figure 14: main figure: Density profiles n(Ep)/A<& in the Abelian 
case a = j6 = 0: for <t>i = 1/157 = 6.410~ 4 (dark red) and for <&2 = 
2/157 = 12.710~ 4 (light cyan). Here the density n(Ep) is divided by 
A<E> = 3>i — S>2> and the plateaus correspond to the energy gaps that 
are highlighted in Fig. 12(b) with the same symbols. Upper right: 
Density profile n(r) in the presence of a typical harmonic potential 
for <J>i = 1/157. The red rectangle specifies the relativistic regime 
which is considered in the main figure. 



Although optical-lattice experiments involving atomic cur- 
rents are on their way, it is preferable to exploit another 
method to evaluate the Hall conductivity. Indeed, several 
works have recently proposed that the quantized Hall con- 
ductivity could be measured through an elegant analysis of 
the density profiles [88, 92-94]. Furthermore, this efficient 
method takes into account the presence of a harmonic poten- 
tial Vharmt/) that generally confines the atoms within the opti- 
cal lattice. If this external potential varies slowly compared to 
the lattice spacing, one may use the local-density approxima- 
tion (LDA) and define a local Fermi energy as 



E F (r) =£>-VharmO)- 



(36) 



The atomic density is then simply obtained through the 
expression n{r) — J deD(e)®(Ep(r) — e), where D(e) is the 
density of states of the uniform system. Each time the Fermi 
energy E-p(r) falls in a gap, the density n(r) depicts a plateau: 
there is a one-to-one correspondence between the spectral 
gaps represented in Fig. 12 and the plateaus of the density 
profiles n(r). We illustrate this phenomenon in Fig. 14, 
where two density profiles for <t>i = 1/157 and O2 = 2/157 
are represented, and where the correspondence with the gaps 
in Fig. 12 (b) has been highlighted with symbols. In this 
figure, we focus on the relativistic region around Ep = and 
set a = j3 = 0. Note that in the following, the density per 
plaquette is studied as a function of the Fermi energy n(Ep), 
being density profiles n(r) simply obtained by considering 
the LDA relation in Eq.(36). In the upper pannel of Fig. 14, 
a typical density profile n(r) obtained for <E>i = 1/157 is 
shown: the relativistic regime can already be identified by an 
unusual behavior (see inside the red rectangle). 



Having identified the unusual plateaus in the density pro- 
files around Ef, one may now evaluate the Hall conductivity 
through the Streda formula [88], which takes the form 



a H (E F ) 



1 dn 
hd$> 



(37) 



First of all, let us show how the Streda relation (37) allows 
one to detect the AQHE at the Graphene point Pi (a = j8 = 
0). Comparing the two density profiles in Fig. 14, the Hall 
conductivity inside the main gaps around Ep = becomes 



Oh (empty rectangle)/; : 
Oh (empty hexagon)/i : 
Oh (empty star)/i : 
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Oh (filled hexagon)/i 
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where ni(£y)/AO [resp. ti2(Ep)/A<f>] corresponds to the 
density profile for Oi = 1/157 (dark red) [resp. O2 = 2/157 
(light cyan)] in Fig. 14, and A<t> = Oi — <I>2. Consequently, 
and in perfect agreement with the Hall plateaus illustrated 
in Fig. 13 (a), the Hall conductivity undergoes the sequence 
Onh £ {...,—10,-6,-2,2,6,10,...}, and therefore agrees 
with the relativistic prediction Oh = ±t(v + 5). Moreover, 
this example shows how a precise analysis of the density 
profile measurements obtained for two close values of the flux 
<t> allows one to observe the AQHE proper to the relativistic 
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Figure 15: Two non-Abelian paths starting from the Graphene point 
(path I) and from the 7T-fiux point (path II). The background contour 
corresponds to |W(cc,j3)| 

regime. Furthermore, these density profiles constitute a pow- 
erful tool to determine the number of Dirac cones situated 
inside the first BZ, since the latter is directly related to the 
jumps observed in the anomalous Hall sequences Nj = 4. 
This important result will be exploited in the following 
section to study how the density profiles are modified as one 
varies the non-Abelian fluxes a and j3. At this point, it is 
worth noticing that the density profiles, and thus the Hall 
plateaus, are not affected as one travels along an Abelian path 
(see Sect. IV B). The latter fact is in agreement with the fact 
that N c / = 4 remains constant along such a path and that no 
topological QPT is thus observed. 



B. The Hall plateaus along a non-Abelian path 

In this section, we consider the two different paths 
within the a — j3 space in Fig. 15. The first path 
starts from the Tl-flux point {a = = tt/2) and reaches 
a = 0.9(371/4), j8 = n/2 in four steps (see the purple dots 
in Fig. 15). The second path starts from the graphene point 
(a = j3 = 0) and reaches a = j8 = n/4 in six steps (see the 
blue dots in Fig. 15). Both paths teach us in what interesting 
manner the density profiles and the Hall plateaus are modified 
as the number of Dirac cones N c / changes. 



Path I.- Let us first study the path starting at the Tl-flux 
point (see path I in Fig. 15). In this path, the parameter 
j3 = k/2 is constant, while a is varied from n/2 to 0.9(37r/4). 
The corresponding energy spectra at vanishing flux <t> = are 
represented in Fig. 16, where one observes how the eight ini- 
tial isotropic Dirac points are progressively deformed as the 
flux is varied. As one travels between the Point 2 (a = 3n/5) 
and the Point 3 ( a = 0.87(3^/4)), one clearly observes 
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(a) (b) 




Figure 16: Energy spectra as one travels along the path I: (a) Point 
1: a = K/2, (b) Point 2: a = 3x/5, (c) Point 3: a = 0.87(3^/4), (d) 
Point 4: a = 0.9(3^/4). The spectra are represented within the first 
BZ and one has set jS = Tt/2 and $ = 0. Dotted circles emphasize 
the merging of two central Dirac points. 

that four Dirac points approach and finally merge at Point 4 
a = 0.9(37r/4) (see Fig. 16 (b) and (c)). Indeed, at criti- 
cality a c = 3n/A, the massless fermions merge into a highly 
anisotropic excitation, the latter being linear in one direction 
and parabolic in the transverse. Below, we explain how this 
topological phase transition is detected through density mea- 
surements and the AQHE. Note that we focus on very low 
fluxes <t>, which allows one to discuss the Hall plateaus in 
terms of the energy spectra in Fig. 17. 

• Point 1: a = n/2. We first compute Oh through the 
density profiles at the Tl-flux point. The van Hove sin- 
gularities (VHS) in the density of states (DOS) delimit 
the relativistic regime between Evhs = ±0.4, where we 
observe steps of eight in the Hall conductivity Oh = 
±|(v+ j). This numerical computation emphasizes 
the presence of the Nj = 8 massless Dirac fermions, as 
observed in Fig. 16 (a). 

• Point 2 to Point 4: a = 3k /5 -> 0.9(3tt/4). As the non- 
Abelian flux a is varied, we observe a transition from 
o H = ±|(v + I) to o H = ±\(v + \) around E = 0. 
Consequently, the density profile measurements per- 
formed along this path witness the topological phase 
transition characterized by the change )V ( / = 8-) 4. 

The above analysis confirms the aforementioned conjecture 
announced in Eq. (33), and thus proves that density measure- 
ments do indeed witness the topological phase transition. 

Path II.- Let us then focus on the second path (see path II 
in Fig. 15). The corresponding energy spectra at vanishing 
flux <t> = are represented in Fig. 17, where one observes that 
the two initial 2-degenerate Dirac points are progressively 
split into four non-degenerate singularities, and at some point 
become anisotropic (see Fig. 17 (d) and (f)). We study below 
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(a) (b) 




Figure 17: Energy spectra as one travels along the path II: (a) Point 
1: a = ft = 0, (b) Point 2: a = P = 7i/16, (c) Point 3: a = ^/8,j8 = 
jr/16, (d) Point 4: a = j5 = jc/S, (e) Point 5: a = k/$,P = x/4, (f) 
Point 6: a = J3 = 71 / 4. The spectra are represented within the first 
BZ and one sets <& = 0. 



how these modifications affect the anomalous QHE. 



• Point 1 , a = /3 = 0: the graphene point has already been 
studied in the previous section, where the anomalous 
steps of four have been observed between EE [—1,1]. 
This effect is due to the presence of two distinct 2- 
degenerate Dirac fermions Nj = 4 at <t> = 0. 

• Point 2, a = [} = ft/ 16: as the degeneracy is lifted, 
four (non-degenerate) Dirac cones are present in the 
Brillouin zone. Our numerical results show that the 
Hall plateaus exhibit steps of four Oh = i^( v + j) 
in the range E G [—0.11,0.11], steps of two in the 
range E € [ — 1, — 0.11)U(0.11, 1], and steps of one 
elsewhere \E\ > 1. These three regimes are delim- 
ited by VHS in the DOS at £ VHS i = ±0.11, and at 
£y HS 2 = ±1. This analysis suggests the existence of 
three different regimes, namely, a relativistic regime 
with four isotropic Dirac cones, an intermediate rela- 
tivistic regime that is characterized by the reminiscent 
of the two initial cones, and finally a non-relativistic 
regime where the linear dispersion vanishes. Addition- 
ally, we have numerically confirmed that the central 
plateau at Ef = (i.e. o# = — 2 — > 2) is robust, which is 
in accordance with symmetry and topological consider- 
ations on the corresponding zero-mode [21]. Therefore, 



the zero-energy Hall conductivity is an ideal witness of 
the number of emerging massless fermions. 

• Point 3 to Point 6, a = = 7r/16-> a = j8 = n/4: 
along this path, one observes the interesting fact that the 
density profiles, and thus the Hall plateaus, are affected 
by the anisotropic character of the Dirac cones. Hence, 
unusual behavior is observed in the relativistic regime 
where the conductivity evolves with steps of two in the 
anisotropic regime (cf. Figs. 17(d) and (f)), whereas 
it evolves with steps of four in the isotropic regime 
(cf. Figs. 17(c) and (e)). One may conclude that the 
anisotropic cones play a peculiar role in the anomalous 
quantum Hall effect: these anisotropic singularities do 
not contribute to the Hall conductivity for Ef ^ 0. 

VI. EXPERIMENTAL CHARACTERIZATION OF 
TOPOLOGICAL PHASES 

In the previous section, we have described in detail a 
method of detection of the phenomena predicted in this paper 
using density profile measurements. This is the easiest and 
the most convenient method, and with current experimental 
setups can be realized with the accuracy of, say, 5 percent. It 
is reasonable to expect that by optimizing this method, accura- 
cies of 1 percent are achievable. However, since the detection 
in quantum many-body systems is one of the most challenging 
tasks, we discuss in this section other possible measurements 
strategies to detect topological phases and Dirac fermions. We 
present the state-of-art of all methods and discuss their preci- 
sion. All the detection methods discussed below have at least 
some of the following properties: 

• Spatial resolution. That means that the methods allow 
to resolve spatial (or equivalently Fourier transformed 
momentum) characteristics of the system. 

• Color resolution. Note that the internal degrees of free- 
dom used for decoding the color states will correspond 
to internal spin, or pseudo-spin states. Many methods 
will allow to resolve color states spectroscopically. 

• N on- demolition property. That means that the methods 
will affect the measured quantum mechanical systems 
in a "minimal fashion", and can be even repeated to al- 
low for time resolution. Only few methods have this 
property, such as contrast imaging of atomic density 
and spin polarization spectroscopy (see below). 

Detection methods can be divided into two groups: those 
which address directly the properties of Dirac fermions, and 
those that have more general character and provide general 
properties of the system. 

a. Direct detection of Dirac fermions These methods 
can be further classified as: i) global detection mechanisms, 
that allow us to infer the total number of massless fermions 
(this is the reason for calling them global), and ii) local detec- 
tion mechanism, that should allow us to map the whole Fermi 
surface, and thus locate each of the conical singularities. 
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Among the global detection mechanisms we can list: 

• Measurement of the atom density as a function of the 
chemical potential: Following [9], a measurement of 
the atom density close to zero chemical potential de- 
pends on the number of Dirac points in the Fermi sur- 
face. In Ref. [9] it was proposed to detect massless- 
to-massive Dirac particles transition; here we propose 
to use it to detect the change of number of the Dirac 
fermions. In principle, this method requires nothing but 
a precise measurement of the density, as in the case of 
the method discussed in the previous section. The semi- 
nal steps toward the precise density measurements were 
achieved quite recently in the first observations of the 
wedding cake structure of the Mott insulator -superfluid 
transition in the optical lattice loosely confined in an 
harmonic trap [95, 96]. Perhaps the most recent state- 
of-art paper is Ref. [97] in which validation of the quan- 
tum simulator of MI-SF transition is realized. Here, 
the most accurate Quantum Monte Carlo simulations 
are directly compared with experimental data; in this 
way accuracies of 1 percent in density determination 
are reached. 

The local detection mechanisms are obviously more de- 
manding; among them we can point out: 

• Atomic angle-resolved photoemission spectroscopy 
ARPES method: a possible candidate for the adapta- 
tion of ARPES techniques to optical lattice setups, is 
based on using Raman transitions to a side level that is 
decoupled from the system [98]. This is a very pow- 
erful method that allows to determine directly the spa- 
tially resolved single particle density matrix, i.e. indi- 
rectly to map the Fermi surface. This new measurement 
technique has been realized for ultracold atom gases in 
Ref. [99]; like photoemission spectroscopy for elec- 
tronic materials, it directly probes low energy excita- 
tions and thus can reveal excitation gaps and/or pseudo- 
gaps with accuracy less than 5 percent. Atomic ARPES 
is a destructive method, but in principle allows for color 
resolution (by splitting the relevant states in an external 
fields and using the spectral resolution), i.e. it allows 
to measure p st (x,x) — ( y ¥ s J \'(x) l l' t (x)), where s,t are the 
color indices. 

• Direct band mapping is an alternative approach [100- 
102]. In this method one releases adiabatically the opti- 
cal potentials to measure directly the quasi-momentum 
distribution. This can be combined with contrast imag- 
ing, so that can be combined with other measurements, 
such as density measurement in the same experimental 
run (see for example [6]). This method allows to map 
higher bands independently, which in principle would 
determine positions of the Dirac points directly. Also, 
it can incorporate the color resolution. 

b. Indirect detection of Dirac fermions Apart from 
these direct methods, there is a whole variety of general de- 
tection methods that will provide additional information and 



characteristics of the system. Although such measurements 
will not give information about the number of Dirac points, 
they will characterize additional interesting aspects of the sys- 
tem: 

• Density imaging: Of course, the density measurements 
are the necessary ingredient of the direct methods based 
on Streda formula, or dependence on chemical poten- 
tial, discussed above. Here, we consider them in gen- 
eral sense as a source of additional information. The 
density imaging allows for color resolution, and can be 
done destructively using absorptive imaging, or non- 
destructively using contrast imaging. It is used in the 
band mapping technique, but can be also used if we turn 
the optical potentials suddenly off to measure the mo- 
mentum distribution, which can give a particularly in- 
teresting characterization of the considered states of the 
system. 

• Noise interferometry: this very powerful (although de- 
structive) method, proposed in Ref. [103], can be used 
to measure directly the two-point density-density cor- 
relation functions, which can be expressed in terms of 
Fourier transforms in the (quasi) momentum space of 
occupations of the corresponding single particle states. 
Originally applied to bosons [104], it is very useful to 
characterize fermionic states [105]. It gives indirect in- 
formation about density of states close to the Fermi en- 
ergy. 

• Hall conductivity measurements: this is discussed in 
detail in Sec. V, where we provide very detailed exam- 
ples for predicted results of such measurements. This 
method is in the first place an indirect method. Still, 
since the anomalous quantum Hall effect is highly af- 
fected by the Fermi surface topology at Ef = 0, the Hall 
conductivity, determined from density profiles mea- 
surements as in Refs. [88, 92-94], provides direct in- 
formation about the change in the number of Dirac 
fermions. 

• Spin polarization spectroscopy: This method (see [106] 
and references therein) may be applied to measure 
Fourier component of the spin density and their fluc- 
tuations. The methods is based on the quantum Fara- 
day effect, i.e. coupling the polarized light beam to 
atoms, mapping the state of atoms onto light, and then 
"reading" off the information about atoms by measur- 
ing the fluctuations of the light beam using homodyne 
detection. It combines color and spatial resolution with 
non-demolition character, and is thus very promising 
and powerful. This method is very well developed ex- 
perimentally (for a review see [107]), and the experi- 
ments with spatial, or, better to say, momentum resolu- 
tion have been initiated in the group of E. Polzik. Sim- 
ilarly, as noise interferometry, it gives indirect informa- 
tion about density of states close to the Fermi energy. 

• Momentum-resolved Bragg spectroscopy: this method 
was pioneered by W.D. Phillips [108] and W. Ketterle 
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[109, 110] groups. It allows to measure the structure 
factor, spin structure factor and other momentum and 
frequency resolved correlation functions. It also gives 
direct information of the dispersion relation of low- 
energy (quasi) particle excitations. The current state- 
of-art of this method is presented in Ref. [Ill], where 
the accuracies of few percents are reached. 

• Methods involving cavity QED: these ideas were 
mostly developed by H. Ritsch and co-workers ([112, 
113], see also [1 14]), but have not yet been realized ex- 
perimentally. They are somewhat similar to ARPES, 
with cavity-enhancing of the signal. In principle they 
are aimed on measuring the spatial density-density cor- 
relations by coupling the systems to the cavity mode, 
and measuring statistics of the light outgoing from the 
cavity. 

Last, but not least, we should mention enormous progress in 
detecting atomic density in situ and in expanding clouds of a 
single atom resolution [115-117]. These methods will add 
on unprecedented precision to many of the above discussed 
detection schemes. 



VII. CONCLUSIONS 

In this article we have discussed the optical-lattice analogue 
of graphene subjected to additional non-Abelian gauge fields. 
We have shown how a long-wavelength perspective captures 
the relativistic nature of the excitations responsible for trans- 
port properties of the system. Besides, when the effective 
fields are switched on, the topology of the Fermi surface is 
modified, and novel phases with gapless fermionic excitations 
that differ from graphene arise. Such topological phases have 
been fully characterized by a topological charge assigned to 
each of the massless fermions. 

We have also discussed how an adiabatic variation of the 
external fluxes along a non-Abelian path can be interpreted 
as a topological phase transition between different quantum 
orders. We have identified two processes responsible for such 
phase transitions, namely, spontaneous and stimulated scatter- 
ing events where the number of gapless fermions is modified. 

We have demonstrated that the creation and annihilation 
processes not only modify the Fermi surface's topology but 
also lead to interesting generalizations of the anomalous QHE. 
Moreover, we proved that a measure of this effect through 



atomic density profiles gives a direct signature of the num- 
ber of Dirac fermions. In this sense, we emphasized that the 
diverse topological phases are characterized by topological 
charges (winding numbers), but also by a specific sequence 
of the underlying Chem numbers. 

These systems provide a very challenging, but at the same 
time rich and interesting playground for applying various 
known detection techniques. Some of them provide direct in- 
formation about Dirac fermions and topological phase transi- 
tion, the other provide additional, but highly relevant charac- 
teristics. Despite the fact that quite a lot of methods exist, it 
is desirable to search for new ones especially suited for the 
systems considered here (attempting for instance to measure 
topological invariants directly). 

We discussed here the physics of single particles. Obvi- 
ously, it opens an avenue towards studies of even more chal- 
lenging physics of interacting ultracold atoms in non-Abelian 
gauge fields, such as the fractional QHE or the Fermi color- 
superconductivity. 

Summarizing, the effect of non-Abelian gauge fields on the 
emerging relativistic theory is far from trivial. Not only do 
such fields transport the topological charges attached to these 
Dirac points, but they also lead to scattering processes where 
pairs of fermions with opposite topological charges are cre- 
ated or annihilated, as a result of the transition between dif- 
ferent quantum orders. The possibility to observe such ex- 
otic phenomena, and extensions thereof, in an optical lattice 
table-top experiment seems to be a challenging and reward- 
ing task where interdisciplinary effects of condensed-matter, 
high-energy, and atomic physics can be observed. 
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